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that time, literally hundreds of crystal 
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using single crystal and powder neutron 
diffraction data, collected with a variety 
of instruments. This work has been 
usually done in collaboration with other 
NBS/NIST divisions and/or universities 
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the technical developments and the 
experimental work, also theoretical aspects 
of crystal geometry have been clarified, 
and significant improvements in the 
techniques of profile refinements have 
been made. It is therefore understandable 
that a comprehensive description of all 
the crystallographic studies carried out up 



to the present is impossible under the 
constraints of space and time imposed by a 
review of this type, and, in the following 
sections, we will limit ourselves to give, 
only a brief account of the topics which, 
in our opinion, represent the highlights of 
the work carried out at the reactor. 
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1. Introduction 



Neutron diffraction at NBS/NIST started soon after 
the NBS reactor became operational in the summer 
of 1969. Since that time, literally hundreds of crystal 
structures have been determined and refined using 
single crystal and powder neutron diffraction data, 
collected with a variety of instruments. This work has 
been usually done in collaboration with other NBS/ 
NIST divisions and/or universities and industrial labora- 
tories. In parallel with the technical developments and 
the experimental work, also theoretical aspects of 
crystal geometry have been clarified, and significant 
improvements in the techniques of profile refinements 
have been made. It is therefore understandable that a 
comprehensive description of all the crystallographic 



studies carried out up to the present is impossible under 
the constraints of space and time imposed by a review of 
this type, and, in the following sections, we will limit 
ourselves to give, only a brief account of the topics 
which, in our opinion, represent the highlights of the 
work carried out at the reactor. 

It is well known that crystal structure analysis 
provides information, which is fundamental for the 
understanding of the physical and technological proper- 
ties of materials. To mention a few examples, drug 
design would be impossible without the detailed 
structural knowledge of the complex molecules involved 
in biological processes, and the assessment of safety 
margins in metallic parts and components does require 
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measurements of the residual stress done with neutron 
and x-ray diffraction methods. As described in other 
chapters of this volume, diffraction measurements are 
also essential in the study of the magnetic behavior of 
many materials, and in clarifying how systems such as 
ice change under the effect of high pressure. 

The instrumentation available to the crystallographer 
today at the NIST reactor must be considered as the best 
in the United States, and, in particular, the new high- 
resolution neutron powder diffractometer is a world- 
class instrument whose usefulness will last for many 
years to come. The crystallographic equipment, coupled 
with other techniques also available at the reactor, give 
to the intelligent experimenter the possibility to deter- 
mine the crystal structure of many materials, even when 
they can be prepared only in small quantities, and to 
relate the structural results to those obtained with other 
techniques. 



2. Instrumentation 
2.1 Introduction 

By inserting a collimator into the face of a reactor, we 
can extract a beam essentially consisting of neutrons 
which have been slowed down by a large number of 
collisions with the moderator and which are in thermal 
equilibrium at the reactor's temperature. The distribu- 
tion of wavelengths in the coUimated beam may be 
obtained by plotting the function 
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Fig. 1. Plot of the function 4>{\)l const, vs A for the reactor equi- 
librium temperature of 350 K. The product (^(A)dA gives the number 
of thermal neutrons of wavelength between A and A + dA emerging 
from an in-pile collimator in one second. 



will diffract, in the best possible geometric conditions, 
a beam of wavelengths comprised in a small interval of 
A centered at 



\ = 2drn sin 0„ 



(3) 



(^(A) = {const. /X')exp[-h^/2mkT\^)] 



(1) 



as function of A [1]. In this expression, h is Plank's 
constant (6.625 X 10"^^ erg s), k is Boltzmann's 
constant (1.381X10"' erg/K), m the neutron mass 
(1.675 X 10"^'' g) and K the absolute temperature. The 
function 0(A) can be defined by saying that <^(A)dA is 
the number of neutrons with wavelengths between A 
and A -I- dA emerging from the collimator in one second, 
and its maximum value occurs for the wavelength 



X = hl{5mkTXy 



(2) 



A plot of </>(A) versus A is shown in Fig. 1 for 
T = 350 K, the equilibrium temperature of our reactor. 
For crystallographic experiments, we must select a 
narrow band of wavelengths from the beam emerging 
from the in-pile collimator, and this task is usually 
accomplished by using as monochromator a single 
crystal of high reflectivity, cut and oriented so that it 



where d^ is the interplanar spacing of the set of planes 
involved in the diffraction process, and IB^ is the take- 
off angle of the monochromator (Fig. 2). The 
monochromatic beam, after passing through a second 
collimator, impinges on the sample, is diffracted in a 
direction which depends on the geometry of the crystal 
and the technique employed in the experiment, and after 
passing through a third collimator, is measured by an 
appropriate detector. The diffracted intensities depend, 
both for the single crystal and the powder methods, on 
the structure and the physical state of the sample and on 
the luminosity (i.e., transmission) of the diffractometer 
used in the experiment. On the other hand, the shapes of 
the observed peaks are given by the convolution of the 
instrumental profile, due exclusively to the geometrical 
characteristics of the diffractometer, and the pure pro- 
file, i.e., the profile that we would observe if the effects 
of the instrument were negligible. In order to under- 
stand the evolution in the design of the diffractometers 
used at different times at the NBS/NIST reactor, we will 
briefly discuss the above concepts and their influence 
on the quality of the collected data. 
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Fig. 2. Schematic representation of the first powder diffractometer installed at the 
BT-1 beam port of the NBS reactor. In its first version, the instrument was 
equipped with just one counter, which was later replaced by a bank of five detectors 
with an angular separation of about 20°. In the figure, C|, C2 and C3 represent the 
in-pile, monochromatic beam and diffracted beam collimators, having angular 
divergences Uu Uj and a^, respectively. The monochromator M is a single crystal 
of copper cut with the reflecting surface parallel to the (110) crystallographic 
planes and oriented to produce the 220 reflection with a take-off angle 2 ft„ of about 
75°. With this configuration the monochromatic beam has a wavelength A of about 
1.5 A. The positions P and A indicate the parallel and antiparallel positions, 
respectively (see text). 



2.2 Instrumental Resolution and Luminosity 

Let us assume that the three coUimators of a single 
crystal or powder diffractometer have horizontal angular 
divergences ai, 0:2 and a^ (these are defined by the 
spacing and the length of the SoUer slits), and that their 
transmission function in the horizontal plane of the 
instrument (which is coincident with the plane of Fig. 2) 
is Gaussian. Let us also assume that the monochromator 
is a mosaic crystal with blocks having a Gaussian spread 
with width at half maximum jSm, and oriented with a 
take-off angle 20^. It has been shown [2-4] that, under 
these assumptions, any diffraction peak centered at the 
angle 2d (where B is the Bragg angle of the peak) has a 
Gaussian instrumental shape with full width at half 
maximum H1/2 given by the equation. 

//i,2 = t/o(tan"0)/tan"0„)±Vo(tan0)/tan0„)H-W (4) 

where Uo, Va and W are function of the angular diver- 
gences «!, 0:2 and a^ and of ;S,„ (in the case of single 



crystal instruments, also the mosaic spread ^c of the 
sample must be taken into consideration). Equation (4) 
is usually written 



Hm = t/tan"0± Viane+W 



where 



U=Uo /tan" 0„, and V=Vq /tan 0„ 



(5) 



(6) 



In Eqs. (4) and (5) the plus and minus signs apply to 
the anti-parallel and parallel positions, respectively, 
(see Fig. 2) and, since V and Vo are always positive, the 
value of H\i2 for the anti-parallel geometry is always 
greater than that obtained when the peaks are scanned 
with the parallel configuration. The differences of the 
full widths at half maximum for the two geometries are 
significant at all values of 6 (and have been fully 
discussed in Refs. [2-4]), and for this reason all two- 
and three-axes spectrometers are designed to work in 
the parallel position. 
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Equation (4) has a minimum at a value of given by 

tan0„,„ = (yo/2C/o)tan0„ (7) 

with a width 

{Hy2f = W- {Vo'n Uo). (8) 

These expressions show that the resolution function of a 
diffractometer can be controlled, within limits, by 
choosing appropriately the horizontal divergences of the 
collimators and the take-off angle of the monochroma- 
tor. This may not be too important in many cases of 
single crystal work, where reflections are rarely close 
enough to one another to cause overlapping, but is 
essential in the case of the powder method, where 
serious overcrowding of reflections occurs even for 
moderately large unit cells. 

2.3 Powder Diffractometers 

The equations discussed before help to explain the 
evolution in the design of the powder diffractometers at 
the BT-1 beam port at our reactor. In the initial 
machine, the monochromatic beam was produced by the 



220 reflection of a flat copper monochromator with 
mosaic spread /3m of about 15 ', oriented with a take-off 
angle 26^ of about 70°. With this configuration the 
wavelength of the monochromatic beam is about 1.5 A, 
i.e., comparable to the wavelength CuK„, extensively 
used in x-ray diffraction experiments. The horizontal 
divergences of the three Soller collimators were of 
the order of 20 ' and the diffracted intensities were 
measured using a single detector in the interval 
5° < 20 < 100°, in steps of 0.05°. 

The advent of the Rietveld method [5,6] in the analy- 
sis of crystal structures had the effect of increasing the 
complexity of the materials that could be studied with 
the powder method, thus requiring resolutions higher 
than those obtainable with our original instrument. This 
was accomplished in the mid seventies by reducing the 
horizontal divergences of the first and third collimator to 
about 10 ', without changing the basic design of the 
diffractometer and the take-off angle of the monochro- 
mator. The loss of luminosity caused by a coUimating 
system tighter than the one in place before, was more 
than compensated by replacing the single counter with 
a bank of five detectors with an angular separation of 
20° (Fig. 2). Examples of the resolution obtainable with 
this instrument are illustrated in Fig. 3 and sh aw that the 
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Fig. 3. Plot of the full width at half maximum (FWHM) vs2d for three configurations of the in-pile, monochro- 
matic beam and diffracted beam collimators and for the 220 reflection of a copper monochromator of mosaic 
spread /3„, = 15 ' and take-off angle 2(?„=75°. The minima of the three resolution functions range from 40° to 60°. 
The lumininosity of the diffractometer for the 10-10-10 and 6-10-6 configurations is 56 % and 21 % that of the 
10-20-10 configuration, respectively. 
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minima of the resolution functions (i.e., the values of 2 
corresponding to the greatest resolution) are all com- 
prised in a region of 20 around 45° to 55°. This is a 
significant disadvantage because, depending on the unit 
cell volume and the chemical and physical nature of the 
sample, severe overcrowding, and consequent over- 
lapping, of the diffraction lines may occur below or 
above this region of 20. In these cases it is clearly 
desirable to have the best instrumental resolution in the 
angular intervals most affected by overlapping, and 
Eq. (7) shows that this can be best accomplished by 
varying the take-off angle of the monochromator. Such 
change, however, must occur without a significant 
variation of the wavelength of the monochromatic beam, 
and with the lowest possible loss of luminosity. 

The instrument best satisfying these often conflicting 
requirements was built following a design proposed by 



E. Prince, of the Reactor Radiation Division, and is 
schematically illustrated in Fig. 4. Its features have been 
described in detail [7] and will only be summarized 
here. With this instrument, the experimenter can mea- 
sure a powder pattern with any one of three monochro- 
mators having take-off angles 2 0^ equal to 75°, 90° or 
120°. The horizontal divergences of the in-pile collima- 
tor can be set at 15 ' or 7 ', while the divergences 02 are 
fixed at the values indicated in Fig. 4 and those of the 
collimators in front of each of the 32 detectors located 
on the recording arm of the instrument are all equal to 
7 '. With these conditions, and assuming the mosaic 
spreads of the monochromators indicated in the figure, 
one obtains the resolution functions illustrated in Fig. 5. 
As expected, the minima of these functions occur in 
different 20 regions. Thus, for the beam generated by 
the 311 reflection of Ge the minimum is located at 
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Cu(311) 
A.=1.540A 
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?L=1.590A 



Range 0°< 20 < 165° 



Fig. 4. Geometrical features of the new, high-resolution, 32-detector, powder diffractometer now operating at the BT- 1 beam port of the NIST 
Center for Neutron Research reactor. This machine is equipped with three focusing monochromators made of iso-oriented sections, as indicated 
in the figure, and with take-off angles 20„ equal to 75°, 90° and 120°. The horizontal divergences Ui and 02 of the first and second collimators 
are indicated in the figure, and the divergence 03 of the third collimator, placed in front of each detector, is 7 '. Since the angular separation of 
successive detectors is 5°, the diffractometer covers a range 0°<2(?£ 160° with an excursion of only 5°. 
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2d ~ 75° and is mostly used for studying materials with 
large unit cells and for magnetic materials, since in 
these cases the intense peaks needed to determine and/ 
or refine the structure are concentrated in the low IB 
region of the pattern. On the other extreme, the 531 
reflection of the Si monochromator, with its minimum 
at about 120°, offers high resolution in the region 
90° < 20 < 130°, and is most effectively used in those 
cases in which it is important to detect splits of the 
diffraction lines (as in the case of phase transitions), or 
when it is essential to measure with high precision the 
thermal factors of the atoms in the structure. Finally, the 
311 reflection of copper is a compromise between the 
two configurations described before, and is routinely 
used to refine structures of average unit cell volumes 
(between 1000 A^ and 4000 A^), such as those of most 
ceramic materials, superconductors, etc. The luminosity 
of the diffractometer is significantly enhanced by the 
use of curved, focusing monochromators and, thanks to 
the presence of 32 detectors, the time required to record 
an entire diffraction pattern, with samples as large as a 
few grams of material, is reduced to a few hours. 

As we have mentioned previously, the observed peak 
shapes are the convolution of the instrumental profile 
and the profile generated by the sample. If both these 
functions are Gaussian with full width at half maximum 
Hii2 and hm respectively, then their convolution is also 
Gaussian, with width Bm given by 



-Bi/2 - H\j2+hm 



(9) 



If hm is small compared to H1/2 (i.e., the size of the 
crystallites in the sample is in the range of 
1500 A to 3000 A with no significant structural distor- 
tions and/or strains), then B^a ~ Hm and Eq. (5) is a 
valid description ofBm as function of 6. If, on the other 
hand, the pure profile is broadened by crystallite size 
effects, then we may express h^a by means of Scherrer's 
Eq. (8) and introduce it in Eq. (5) obtaining the expres- 



B\i2 = [U+(K^\yD^)tiin^e]-Vtane+[W+(K^\VD^)] 

(10) 

where ^ is a constant (—0.9) and D is defined as the 
thickness of the crystallites in the direction perpendicu- 
lar to the reflecting planes (it should be noted that, sensu 
s trie to, Eq. (10) is valid only when the instrumental and 
pure profiles are both Gaussian). When the crystallites 
have odd shapes (such as plate-like or needle-like), D 
changes from reflection to reflection and an expression 
like Eq. (5) cannot be used to describe the dependence 
of Bi/2 on 0. If, on the other hand, the crystallites are 
approximately spherical and of uniform size, D is con- 
stant and Eq. (10) has the same form as Eq. (5). It is 
remarkable that in the majority of cases encountered in 
practice, this condition is satisfied and the observed 
peaks are very nearly Gaussian with widths having a 
dependence on well described by Eq. (5). Examples of 
single peaks obtained with the high-resolution diffrac- 
tometer now in use at the BT- 1 beam port are shown in 




^<' 120° (Si 531) 



Fig. 5. Examples of resolution functions for three monochromators of our new, 
high-resolution diffractometer. The take-off angle 2(?„, and the monochromator 
characteristics are indicated in the figure. It should be noted that, consistently with 
Eq. (7), the minima of these functions depend on the value of S„ and that, in this 
case, they occur at values of 20 about equal to the values of 20,„ of each function. 
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Fig. 6. The excellent quality of the data that can be 
collected with the instrument is clearly illustrated in 
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Fig. 6. Gaussian fit of the 222 reflection obtained from a powder of 
the cubic perovskite superconductor MgCNis . The crosses are exper- 
imental values and the continuous line the calculated profile, and the 
difference between the two is shown at the bottom of figure. The good 
agreement between observed and calculated intensities is proof of the 
high quality of both the new diffractometer and the sample used in the 
experiment. 



Fig. 7, which shows the agreement between observed 
and calculated intensities after the refinement of the 
structure of MgCNi3 [9]. It is proper here to point out 
that this unusual and important perovskite compound is 
a superconductor with Tc ~ 9 K, and that its structure 
was determined and refined by Q. Huang, a guest 
scientist of the NIST Center for Neutron Research, using 
a sample prepared in the Department of Chemistry of 
Princeton University. 

2.4 Single Crystal Diffractometers 

Since overlapping of reflections is rarely a problem in 
single crystal diffraction work, the neutron diffracto- 
meters used in this type of analysis are usually designed 
to maximize luminosity and without paying too much 
attention to resolution. Consistently with this concept, 
the single crystal instruments installed at the BT-7 and 
BT-8 beam ports used graphite or copper monochroma- 
tors with large mosaic spreads and essentially no colli- 
mators. The crystal mounting is in general a three-circle 
goniometer of the type schematically represented in 
Fig. 8. The rotations d, x ^iid "^ allow the experimenter 
to bring any reflection h k I into diffracting position on 
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Fig. 7. Calculated (continuous line) and observed (crosses) patterns, plotted after the final refinement of the cubic perovskite 
superconductor MgCNis . The difference between theoretical and experimental values is indicated by the profile in the lower 
part of figure. The short vertical marks indicate the Bragg positions of the reflections of a small quantity of graphite impurity 
(high marks) and of the MgCNii sample (lower marks). 
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Fig. 8. Goniometer head of a four-circle diffractometer used for single crystal diffraction experiments. In this 
instrument the detector rotates about the 20 circle axis, and the crystal is located the center of the x circle. 
Intensities are measured on the horizontal plane, which is parallel to the planes of the fl and 20 circles. The 
ft) rotation (which takes place about the 20, but is independent from the $-20 rotation) allows the experi- 
menter to rotate the crystal about the scattering vector of any of the measurable reflections. 



the equatorial plane, defined by the incident and 
diffracted beams, and to rotate the crystal of any angle 
\\i about the scattering vector. This last option is impor- 
tant in those cases in which multiple diffraction must be 
avoided [10]. Our original diffractometer was equipped 
with a single detector rotating about the axis of the 
20 circle, and the diffracted intensities were mea- 
sured following the automatic procedure described in 
Ref. [11]. 

During the seventies, the BT-7 beam port was used to 
develop a new type of single crystal instrument suitable 



for protein crystallography. Since the unit cells of 
protein materials are very large, it is essential to design 
the instrument so that as many reflections as possible 
are measured simultaneously. Consistently with this 
requirement, our diffractometer was equipped with a 
linear position- sensitive detector, and the crystal, 
mounted on a three-circle goniometer of the type 
illustrated in Fig. 8, was rotated according to the 
flat-cone Weissenberg geometry [12]. With this 
technique, all reflections of any given layer of the 
reciprocal lattice develop on a single plane, which is 
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made to coincide with the plane formed by the linear 
detector and the incident beam, as shown schematically 
in Fig. 9 [13]. This machine became obsolete with the 
advent of two-dimensional detectors, but in its lifetime 
valuable data were collected with it, which were used to 
solve the structure of the bovine pancreatic trypsin in- 
hibitor [14]. 
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Fig. 9. Geometry of the flat-cone Weissenberg method. The c-axis 
(or any other important rational direction of the crystal) is first 
brought into the plane of figure, at the position indicated, by a go- 
niometer of the type illustrated in Fig. 8. In this position the l-th level 
of the reciprocal lattice, which is perpendicular to the c -axis, passes 
through the center of the reflection sphere. When the crystal is rotated 
about the c-axis by using the appropriate combinations of the <^, x 
and to rotations, the reflections of the / -th level develop on a vertical 
plane (also perpendicular to c) and can be intercepted by a linear 
position sensitive detector (whose projection on the plane formed by 
the primary beam and the c-axis is shown in the figure). 



3. Structural Studies 
3.1 Neutrons and X-Rays 

In order to understand the type of structural problems 
studied by neutron diffraction at the NBS/NIST reactor 
is it necessary to briefly summarize the main properties 
of the neutron radiation and to underline the differences 
between the neutron techniques and the older and more 
conventional x-ray methods. 

The wavelength of monochromatic beams of thermal 
neutrons used in structural investigations is of the order 
of 1 A, i.e., similar to that of x rays used for the same 
purpose. The two types of radiation, however, are scat- 
tered by atoms in quite different ways. In the case of 



X rays, the fundamental scattering is generated by 
the extra nuclear electrons, which, by virtue of their 
charge, interact with the incident beam. In the case 
of neutrons,the fundamental scattering is generated by 
the atomic nuclei, except for magnetic materials where 
also electronic scattering can be significant (see 
article entitled "Magnetic Structure Determinations at 
NBS/NIST" in this Special Issue). The two scattering 
processes result in qualitative and quantitative differ- 
ences in the way in which structural experiments are 
carried out and in the type of results obtained in the two 
cases, as we will soon describe. 

Since the dimension of the electron clouds surround- 
ing the atomic nuclei are comparable with the wave- 
lengths of the x-rays used in crystallography (0.7 A 
to 2.0 A), the amplitude of the scattered radiation will 
fall off rapidly when the angle between incident and 
scattered directions increases. On the other hand, 
neutron scattering is essentially isotropic because the 
size of the atomic nuclei is small compared with the 
wavelength of thermal neutrons. In other words, the 
scattering amplitude of x rays is function of the diffrac- 
tion angle, while that of neutrons remains practically 
constant with 0. This behavior is illustrated in Fig. 10, 
where the x-ray and neutron scattering amplitudes of an 
yttrium atom are shown as function of (sin0/A). The 
consequence of this effect is that the intensities of x-ray 
diffraction patterns are significantly weakened at high 
values of 6, while those of neutrons are not. In practical 
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Fig. 10. Plot of the neutron and x-ray scattering amplitudes of 
yttrium versus (sin 0)/A. 
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terms this means that in a neutron experiment we can 
collect data also at high values of , where the intensi- 
ties are very sensitive to both the positional and thermal 
parameters of the atoms in a structure. It is worthwhile 
to note that a precise knowledge of the atomic tempera- 
ture factors is important not only to evaluate appropriate 
corrections to the bond lengths, but also to understand 
the role of molecular motions in phase transitions. An 
example of this application of neutron diffraction is the 
study by Choi and Prince of several metal azides 
(specifically, fi-NaNs, KN3, RbNj, and TIN3) done at 
the NBS reactor [15,16]. In this series of experiments it 
was shown that the azide ions are linear and centro- 
symmetric and that their thermal vibrations could be 
treated as rigid body motions. 

As mentioned previously, since the x rays interact 
with the extra nuclear electrons, the amplitude scattered 
by an atom is generated by the contribution of all the 
electrons surrounding the nucleus. It follows that the 
intensities diffracted by crystals containing heavy and 
light atoms at the same time will be dominated by the 
heavy atoms. While this property has been used in the 
past as a way to solve a large number of crystal struc- 
tures (especially, but not exclusively, of organic com- 
pounds) [17-19], in many cases its presence makes it 
difficult to locate with precision the light atoms and to 
determine their thermal and occupancy parameters. In 
the case of neutrons the scattering amplitude of all 
atoms is in general of the same order (within a factor of 
2 or 3) and does not show a regular dependence on the 
atomic number. As a consequence, positional and 
thermal parameters of light atoms can be determined in 
a neutron experiment even when the number of heavy 
atoms in a compound is large. Examples illustrating this 
property of neutrons are numerous in the literature. In 
the case of the azides mentioned previously, the position 
of the N atoms could be easily determined in the neutron 
experiment since the coherent scattering length b of 
N is of the same order as those of the metal atoms 
(e.g., 0.926 X 10" and 0.878 X 10'" cm for N and Tl, 
respectively; note that the atomic numbers of the 
two elements are 7 and 81). Another example is the 
123 superconductor YBa2Cu306+.t (with ^x< 7), for 
which the scattering lengths of oxygen, copper, yttrium 
and barium are 0.581 XlO-'^ 0.772X10", 0.775X 
10 ", and 0.507 X 10 " cm, respectively, and the corre- 
sponding atomic numbers 8, 29, 39, and 56. This situa- 
tion has allowed crystallographers to solve the structure 
with great precision and to determine the number and 
the location of the oxygen vacancies over the entire 
range of x [20-22]. 



An additional, important difference between x rays 
and thermal neutrons is the extent to which they are 
absorbed in materials. In a significant number of cases, 
which include the majority of heavy atom elements, the 
absorption of thermal neutron is negligible, while that of 
X rays is very large and requires special experimental 
procedures to reduce its effects and/or the evaluation of 
complicated corrections to the observed intensities, even 
when small samples are used. This point is illustrated 
in Table 1 for the 123 compound YBa2Cu307. The 
absorption of this material for x rays and neutrons 
of wavelengths 1.54 A and 1.08 A, respectively, is 
expressed in terms of the mass absorption coefficient 
(jil p) (where jU, is the linear absorption coefficient in 
cm ' and p density of the sample in gem "*) because this 
quantity is approximately independent of the physical 
state of the compound and its value can be obtained, to 
a good approximation, by adding the mass absorption 
coefficients of the elements according to the equation 



l±lp = Xiy,{iilp)i 



(11) 



where y, is the mass fraction contributed by element / of 
mass absorption coefficient {ill p ),, and the summation 
is taken over all the constituent elements /. The low 
value of /x/p evaluated for neutrons means that the 
attenuation of the primary beam in a neutron experi- 
ment is essentially due to the diffraction taking place in 
the sample (secondary extinction), rather than to true 
absorption. An important consequence of this property 
is that it makes it possible to use large or bulky samples, 
thus reducing the data collection time, and to minimize, 
or eliminate, preferred orientation effects in experi- 
ments done with the powder technique. As we will show 
later, the low absorption of thermal neutron has been 
used to great advantage in texture studies of bulky 
samples of metals and alloys. 



Table 1. Absorption in YBa2Cu307 (molecular weight: 666.5) 



Elements Atomic wt. 



X-Ray 

(A= 1.546) A 

[Lip ii./py' 



Neutron 
(A= 1.08) A 



(j,/p 



iL/py 



Y 


88.9 X 1 


134 


17.9 


0.006 


0.001 


2Ba 


137.4 X 2 


330 


136.1 


0.003 


0.001 


3Cu 


63.6 X 3 


52.9 


15.2 


0.002 


0.001 


70 


16.0 X 7 


11.5 


1.9 


0.000 


0.000 


2 


666.5 




171.0 




0.003 



" cm' g-' [23] 

'' -y is the mass fraction of the corresponding chemical element. 
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3.2 Crystal Structure Analyses 
3.2.1 Solid State Ionics 

In general, the structures of ionic conductors are 
made of two parts, one formed by atoms or ions in fixed 
positions and the other by ions of high mobiUty, which 
account for the ionic conductivity. The atoms in fixed 
positions form the so-called framework structure [24], 
which may have a wide range of complexity. For exam- 
ple, in a-AgI it is a simple body centered arrangement 
of I ions with cubic symmetry [25]; in Na(3- and 
(3"- alumina it is a block of two or three units of a 
spinel-like structure [26]; and in some zeolites it may be 
a still more complex assemblage of cage-like units [27]. 
An accurate knowledge of the geometry of the frame- 
work structure is essential for understanding the con- 
duction mechanisms in these compounds. In fact, the 
mobility of the disordered ions between sites depends, 
among other things, on the size of the opening through 
which the ions move. An example of the analysis of the 
geometrical factors involved in this process is given in 
Ref. [28] for the compound Na4Zr2Si30i2. In this struc- 
ture there are theoretically four possible pathways for 
diffusion, but only one of them is sufficiently large 
to allow reasonable mobility of the Na* ions at room 
temperature. 

The mobile ions are somehow distributed among sites 
of multiplicity greater than the actual number of ions in 
the unit cell, and their diffusion within the frame- 
work structure may occur with different mechanisms. In 
a-AgI, for example, the two Ag* cations in the unit cell 
are distributed over twelve energetically equivalent sites 
and may move in a three-dimensional system of path- 
ways [25]. In sodium /3"-alumina, the Na* cations are 
free to move in two-dimensional layers between the 
rigid spinel-like blocks [26]. Finally, in materials such 
as ZrOz stabilized with CaO, the O^" ions are located 
in well defined crystallographic sites which are only 
partially occupied, and the conduction takes place 
through this defective structure, with the oxygen ions 
diffusing from the occupied to the vacant sites [29]. 

The above discussion shows that, in order to prepare 
compounds which are potential ionic conductors it is 
necessary to start with appropriate framework struc- 
tures into which the mobile ions are inserted. The inser- 
tion reaction changes the chemical and structural nature 
of the host in ways that are largely unpredictable. A 
knowledge of the insertion mechanism and of the 
changes that the inserted ions cause is therefore needed 
to prepare solid state ionics and/or to predict which 
framework structures are most favorable for the forma- 
tion of these useful materials. In particular, compounds 
which undergo topotactic insertion of lithium are inter- 
esting because of their potential use as electrodes in 
secondary batteries. Lithium is ionic in these materials 
and the charge introduced by Li* in the host structure is 



compensated by a corresponding reduction in the oxida- 
tion state of the host cations. Since the insertion of extra 
ions usually causes significant changes which may 
pulverize the starting materials, it is in general difficult 
or impossible to grow single crystals of the inserted 
compounds, and consequently the powder technique, 
coupled with profile refinement analysis, is the only 
method for studying this type of structures. As we have 
discussed previously, neutrons must be used, since they 
are more sensitive than x-rays in locating the light Li* 
ions in the structure. 

Lithium inserted materials were being studied in the 
early 1980s at the Ceramics Division of NBS and at Bell 
Laboratories' in Murray Hill, N. J. It was therefore 
natural for the Reactor Radiation Division to establish 
collaborations with these organizations, and a number of 
problems were studied and solved over the years under 
this project. In what follows we will briefly discuss, as 
an example of the work done, the results obtained for the 
system LivReOj. 

Compounds such as ReOs (and WO3) are suitable 
framework structure for Li* and, in fact the system 
Li^Re03 exists in three distinct phases with composition 
< X < 0.35, X = 1.0 and 1.8 < X < 2.0. The structure 
of ReOs has the perovskite arrangement of the ReOe 
octahedra (Fig. 11) and the cuboctahedral cavity, nor- 
mally occupied by the large A cations in the A B O3 
perovskites, is in this case empty. For concentrations 
less than 0.35 Li atoms per formula unit, the Li* ions 
occupy this cavity [30]. The symmetry for the composi- 
tion Lio.2Re03 is that of the cubic space group Im,, and 
is produced by a distortion of the original ReOs struc- 
ture caused by a cooperative rotation of the ReOe octa- 
hedra, as illustrated in Fig. 12. The nature of the distor- 
tion can be understood if we remember that in the 
perovskite configuration (see Fig. 11) the octahedra 
share corners, so that if any one of them is tilted, all the 
others will also rotate. This situation is illustrated in 
Fig. 13 for the simple case in which the tilt axis is 
parallel to the vertical axis of the structure. In Lio2Re03, 
the geometry of the final configuration is more complex 
because a tilt of the same angle occurs about all three 
crystallographic axes simultaneously. It can be shown 
[31] that the tilt angle (p is related to the a -parameter of 
the Im3 structure {a = 7.3979 A) and to the Re-0 bond 
distance {d= 1.887 A) by the equation 



cos (p = (3fl-4 d)/8d 



(12) 



so that, for our structure the tilt angle is cp ~ 14°. It 
is worthwhile noting that in general, in this type of 



Certain commercial equipment, instruments, or materials are identi- 
fied in this paper to foster understanding. Such identification does not 
imply recommendation or endorsement by the National Institute of 
Standards and Technology, nor does it imply that the materials or 
equipment identified are necessarily the best available for the purpose. 
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{A}0,, {B}0, {0}A,B, 



B 



Oo 



Figure 11. The unit cell of the ideal cubic structure of an ABOi perovskite. The large A cation is at the center 
of a cuboctahedron and the small B cations at the center of regular octahedra. In order to exist, the A-O and B-O 
bond distances must be related so that (A-O) = (B-0)(2)"^. On the right side of figure are indicated the bonding 
scheme and the coordinations of the atoms (atoms in curly brackets are the central, coordinating atoms). 





Fig. 12. The cubic (/m3) structure of Lio.2 ReOj projected on the 
(001) plane. The empty triangles represent Re atoms at z = 1/4, 3/4; 
empty squares: Li at z = 1/2; empty circles: Li at 0, 1; full trian- 
gles: O at z= 1/4, 3/4; full squares: O at z= 1/2; full circles: oatz=0, 
1 . In this structure the ReOf, octahedra are tilted of the same angle, and 
in the same sense, along all three crystallographic axes. 



Fig. 13. Example of cooperative tilting of the BOs octahedra in the 
perovskite structure. In this case the octahedra are rotated only about 
the direction of the vertical axis and in opposite senses in each succes- 
sive layer. When tilting occurs, the twelve-fold coordination of the 
large A cations changes drastically, while the BOs octahedra in many 
cases remain practically undistorted. The oxygen and B atoms are 
located at the corners and the centers of the octahedra, respectively. 
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distortions the octahedra remain essentially the same. 
For example, the Re-0 distances in Re03 and in 
Lio2Re03 are practically unchanged, the difference be- 
ing of the order of one hundredth of an A, and the 
0-Re-O angles are quite close to 90°. 

The tilt system that accompanies the Li insertion 
when LiReOs and Li2Re03 are formed is more complex 
and more extensive than the one discussed previously 
[32]. In this configuration, the angle <p is of the order of 
30° and in the resulting structure the cub-octahedral 
cavity of Re03 is changed into two face-sharing octa- 
hedra with six tetrahedra on the re-entrant faces (see 
Fig. 14). The symmetry of this structure is 7? 3c and the 





two octahedra are fully occupied by Li in Li2Re03 and 
only half occupied in LiRe03. This result shows that the 
tilt system of the Re03 host structure produces a coordi- 
nation favorable to the insertion of one and two atoms of 
Li per formula unit, and it is made possible by the fact 
that the ReOe octahedra share corners. 

3.2.2 Crystallography of Superconductors: 
the System Ba (Pbi.^BiJOj 

It is important to point out that high-temperature su- 
perconductivity in ceramic materials was first observed 
in the system Ba (Pbi_j^Bij^)03 [33]. Soon after this dis- 
covery, structural studies were carried out over the en- 
tire range of compositions ^x< 1 [34-39]. Although 
the general structural features of Ba (Pbi_tBi.v)03 do 
not change dramatically with composition, the type of 
distortions present in the basic atomic arrangement does 
vary with x. At room temperature, the symmetry 
changes according to the sequence 



0<x<0.05 
0.05<x<0.35 
0.35<x<0.90 
0.90<x<1.00 



orthorhombic 
tetragonal 
orthorhombic 
monoclinic 



Fig. 14. Schematic representation of the distortion of the cuboctahe- 
dral cavity (top figure) caused by the tilting of the ReO(, octahedra, 
generating two face-sharing octahedra containing Li atoms in the 
structures of LiReO, and LiiReOs. The oxygen atoms are represented 
by the empty circles. 



Superconductivity exists only in the tetragonal phase. 
The value of the critical temperature T^ increases with 
X, reaches a maximum of Tc ~ 13 K for x ~ 0.25, and 
then decreased between 0.25 and 0.35. The arrangement 
of the atoms in this system is closely related to that of 
perovskite and the structure of the superconducting te- 
tragonal phase, schematically shown in Fig. 15, can be 
generated by tilting consecutive (Pbi_;iBit)06 octahedra 
of the same angle (about 8° forx = 0.25), but in opposite 
directions, about the c -axis of the unit cell. 

For a complete characterization of this system we 
have to know the oxidation state of the Bi atoms, since 
two formulations are possible as indicated below 

Ba (Pbfi Bi**)03 and Ba [Pbfi (Bifc BiE:'5)x]03 . 

Most of the structural studies carried out to determine 
the valence of Bi have been done with samples of 
BaBi03 [37,39], in which Bi is located in two crystallo- 
graphically independent positions of the monoclinic 
space group I2/m. Powder neutron diffraction results 
have shown that the distances Bi(I)-0 and Bi(2)-0 are 
significantly different (2.283 A and 2.126 A, respec- 
tively) and, from these values, it was concluded that the 
site Bi(l) is occupied by Bi^* cations, and site Bi(2) by 
Bi^*. X-ray photoemission spectroscopy on single crys- 
tal of BaBi03confirmed that there are two independent 
positions for Bi [40]. However, the difference in elec- 
trons density was found to be too small to justify a 
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O Pb (0.75%)+ Bi (0.25%) 
% Ba 

Oo 



BaPbo.75Bio.25O 
R.T. 14 /mem 



Fig. 15. Crystal structure of the compound Ba(Pbo.75Bio.25)03. Two unit cells are joined together to show the close 
relationship with the unit cell of perovskite, defined by the axes flp, b^ and Cp. The two octahedra shown in the figure are 
rotated of the same angle and in opposite senses a bout the t-axis, according to the scheme shown in Fig. 13. 



complete disproportionation into 3+ and 5+ valence 
states of Bi. This contradiction was resolved in a joint 
study at the NBS reactor and at Bell Laboratories. In 
this research [41] it was shown that the previous result 
may be explained by assuming a different degree of 
ordering of the Bi"** and Bi''* cations located at the two 
inequivalent positions, in the sense that one position is 
occupied by somewhat more Bi^* than Bi^* and the other 
position is occupied in the opposite way. These results 
were corroborated in subsequent neutron diffraction 
experiments at the reactor [42] which showed that there 
are in fact two polymorphs of BaBiOs, one in which the 
site Bi(l) has a valence of about 3.5+ and site Bi(2) a 
valence of 4.5+, and the other in which the two 
sites have both valence 4.0+, i.e., contain an almost 
equal proportion of Bi^* and Bi^*. The charge ordering 
between the two sites, i.e., the existence of one 
polymorph or the other, depends on the method of 
preparation of the samples and on their thermal history, 
and in fact, the two phases can be transformed into one 
another by heat treatment. It is interesting to note that 
the same assumption of a disordered distribution of the 
Bi ions has been made to explain the behavior of the 
system Ba,_^K,-Bi03 [43,44]. 



3.2.3 The System Laj-. M, Cu04-,^ (M = Ba, Sr) 

The discovery of superconductivity in this system was 
made on a sample consisting of a mixture of phases 
[45], and only in subsequent experiments was the super- 
conducting material identified as La2_tBatCu04_^ 
[46,47]. This compound has the K2NiF4-type structure 
at room temperature, consisting of alternate layers with 
the perovskite and rocksalt structure, as illustrated in 
Fig. 16. The value of the critical temperature T^ is also, 
in this case, a function of x and reaches its maximum of 
about 35 K for x^O.15. The detailed structure for this 
composition was determined by neutron powder diffrac- 
tion [48,49] which showed that the symmetry at room 
temperature is tetragonal lAlmmm and that a transition 
to an orthorhombic phase with space group Cmca oc- 
curs at about 180 K [49]. The lattice parameters of the 
two phases are related so that ao~Co~at Vz and fco~c,, 
where a^, bo, and Co are the parameters of the or- 
thorhombic cell and a, and c, these of the tetragonal cell. 

At this stage attempts were being made to prepare 
compounds with values of T^ higher than 35 K, and the 
system La2-j:Baj^Cu04_v was considered as a candidate 
at Bell Laboratories. The composition x~0.15 was 
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.[(CuO,)J/.0),(/.0)„(CuO,)J/.0)J/.0)J(CuO,)„, 



{/.}0(1)A2). 
{Cu}0(l)A2), 
{0(l)}/.,Cu 
{0(2)}/.,Cu, 



Fig. 16. Layer sequence, bonding scheme and atomic coordinations 
in the structure of (Laa-^^M,) Cu04_v (M = Ba, Sr). The symbol L 
represents the mixed cations Lai_5M5, Where 8 = xl2. In the se- 
quence, the formulas enclosed in the parentheses give the chemical 
composition of each layer and the subscripts o and c indicate if the 
cations are at the origin or in the center, respectively, of the mesh 
characteristic of the layer. The square brackets include the content of 
one unit cell of the structure. 



found to have Tc ~ 40 K and its structure was refined at 
the NBS reactor using neutron powder diffraction 
data [50]. This analysis showed that the Ba and Sr 
compounds are isomorphic at room temperature and 
that also the Sr compound undergoes a phase transition 
from tetragonal to orthorhombic at about 200 K. At all 
temperatures, the lanthanum and strontium atoms were 
found to be distributed at random over the same equiva- 
lent positions. A schematic picture of the tetragonal 
structure is shown in Fig. 17, where the perovskite-like 
planes are emphasized to illustrate the two-dimensional 
nature of the copper-oxygen layers perpendicular to the 
c-axis. These layers are separated by the La/Sr-oxygen 
planes, so that the copper atoms in one plane do not 
share oxygen atoms with copper atoms located on other 
planes. The Cu-0 distances within the perovslcite-like 
planes are short (about 1.89 A) while those in the 
perpendicular direction are long (about 2.41 A), and 
consequently the coordination polyhedron around 
copper is a bi-pyramid, rather than an octohedron. The 
shape of the bi-pyramid changes in a subtle way when 
the symmetry is lowered from tetragonal to orthor- 
hombic. More specifically, the base becomes rectangu- 
lar and the copper-oxygen planes are buckled as a 
consequence of small but significant rotations of the 
bi-pyramids about axes parallel to the a- and b- direc- 
tions, as shown in Fig. 18. The neutron diffraction stud- 
ies of the Ba and Sr compounds showed no oxygen 
deficiency. This means that the charge compensation for 
the substitution of La'* with Ba^* or Sr^* must be accom- 
plished entirely by oxidation of copper from 2+ to 3-I-. 
By assuming a simple ionic model, the oxidation state of 
copper in Lai 85M0 15CUO4 is 2.15-I-, i.e., the Cu sites are 
occupied by Cu^* and Cu'* in the ratio 17/3. 




ot 



Fig. 17. View of the room temperature tetragonal structure of 
Laj 35Sro.i5Cu04. The large shaded circles indicate the mixed cations 
Lao925Sr().o75, the small filled circles represent copper atoms at the 
center of bi-pyramids and the oxygen atoms are located at the comers 
of the polyhedra. 



3.2.4 The System YBajCujO, 

Superconductivity with 7'c='94 K was discovered in a 
sample which was a mixture of phases [51] and, again, 
the superconducting compound was identified later in 
several laboratories as YBa2Cu30i, with x~l .Q [52- 
54]. Profile analysis based on neutron powder data 
quickly gave an accurate description of the structure 
which is shown in Fig. 19, next to the structure of 
perovskite to underline the similarities and the differ- 
ences between the two atomic arrangements. The 
refined structural parameters obtained in some of the 
most accurate neutron diffraction studies of the compo- 
sition x~7.0 are given in Table 2, which has been taken 
from Ref. [55], and the bonding scheme and the atomic 
coordinations are illustrated in Fig. 20. 

Figure 19 shows that the structure of YBa2Cu307 is 
obtained from that of perovskite by tripling the c-axis, 
by eliminating in an orderly fashion some of the oxygen 
atoms, and by ordering the Y and Ba atoms along the 
c-axis of the unit cell. The copper atoms are located on 
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o Cu 
^ 0(up) 

O O(down) 






Fig. 18. Schematic representation of the distortion of the copper-oxygen layers in the orthorhombic phase of 
Lai 85Sro 15CUO4. The central part of the figure is the projection of one CuOi layer on the plane of the Cu atoms 
along the c-axis, and the figures at the bottom and the left sides are projections along a, and b,, respectively. 



two crystallographically independent sites of the or- 
thorhombic space group Pmmm. The first, labeled 
Cu(l), has four-fold planar coordination in which the 
near square CUO2 units share corners and form chains 
along the b -axis, and the second, Cu(2) has a pyramidal 
coordination with two Cu(2)-0(2) and two Cu(2)-0(3) 
short bonds (1.930 A and 1.961 A, respectively) and one 
Cu(2)-0(1) long bond (2.295 A). Because of this pecu- 
liar structural feature, there are in the structure two- 
dimensional layers of composition CUO2 perpendicular 
to the c-axis. The oxygen atoms in the layers are slightly 
shifted from their ideal perovskite positions, producing 
a buckling of the Cu(2)-0(23) bonds, as indicated in 
Fig. 19. 

It was soon observed that the oxygen stoichiometry in 
YBa2Cu30x may be lower than seven atoms per formula 
unit, depending on the preparation and the thermal 
history of the sample, and that the oxygen vacancies 
associated with the decrease of the oxygen content in- 
volve only the 0(4) sites of the CUO2 chains [20]. When 
the number of vacancies becomes sufficiently large, the 
sites 0(5), which are empty for x = 7.0, become gradu- 
ally filled and, at some point, the structure from or 



thorhombic Pmmm becomes tetragonal PAImmm. The 
transition is a function of temperature and oxygen 
content, and, using appropriate procedures, the 
orthorhombic modification can be retained at composi- 
tions as low as X = 6.3. The value of the critical temper- 
ature Tc decreases as x decreases, and it becomes equal 
to zero when the orthorhombic to tetragonal transition 
takes place. When x = 6.0, all 0(4) and 0(5) sites are 
empty, and the coordinations of Cu(l) and Ba became 
two and eight-fold, respectively [58]. 

The structural changes occurring in YBa2Cu30x over 
the entire range of compositions 6.0^x<7.0 have been 
studied in two accurate and extensive sets of neutron 
diffraction experiments [21,22], and the interpretation 
of these observations has been provided by Brown, 
using the bond valence method of analysis [59,60]. By 
comparing the a -axis parameter measured experimen- 
tally with that calculated from the theoretical Ba-O(l), 
Cu(l)-0(4) and Cu(2)-0(23) distances consistent with 
different valences of the chain and pyramidal copper 
atom. Brown was able to prove that the oxidation states 
of Cu(l) and Cu(2) in YBa2Cu306 are \+ and 2+, 
respectively, while they are approximately the same 
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Oi 0(3) 




c 



0(5) 
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Fig. 19. View of the structure of the 94 K superconductor Ba2YCu307 (right side) and of the configuration produced by 
stacking along c three unit cells of the ideal structure of perovskite ABO3 (left side). The atomic arrangement of the 
superconductor differs from that of perovskite because (i) the Y and Ba atoms are ordered; (ii) the Y layers and the sites 
0(5) are not occupied by oxygen atoms; (iii) the CUO2 layers are buckled. In Ba2YCu307, the Cu(l) atoms are in four-fold 
planar coordination, resulting in the characteristic chains which propagate along the b-axis of orthorhombic unit cell. 
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Table 2. Comparison of Rietveld refinements for "yBa2Cu307 obtained in four laboratories 



Atom 



Parameter 



Ref. [55] 



Ref. [56] 



Ref [57] 



Ref. [20] 



Ba 


z 


0.1839(3) 


0.1841(3) 


0.1843(3) 


0.1839(2) 




B 


0.4(1) 


0.6(1) 


0.54(5) 


0.65(5) 


Y 


B 


0.2(1) 


0.6(1) 


0.46(4) 


0.56(4) 


Cu(l) 


B 


0.2(1) 


0.4(1) 


0.50(5) 


0.55(4) 


Cu(2) 


z 


0.3546(2) 


0.3549(3) 


0.3556(1) 


0.3547(1) 




B 


0.3(1) 


0.5(1) 


0.29(4) 


0.49(4) 


0(1) 


z 


0.1589(3) 


0.1581(4) 


0.1584(2) 


0.1581(2) 




B 


0.5(2) 


0.9(1) 


0.67(5) 


0.78(5) 




n 


2.01(3) 


2.0 


2.0 


2.0 


0(2) 


z 


0.3783(3) 


0.3779(4) 


0.3773(2) 


0.3779(2) 




B 


0.5(2) 


0.1(1) 


0.56(5) 


0.57(5) 




n 


2.01(3) 


2.0 


1.89(2) 


2.0 


0(3) 


z 


0.3780(3) 


0.3777(5) 


0.3789(3) 


0.3776(2) 




B 


0.4(2) 


0.3(1) 


0.37(5) 


0.55(5) 




n 


2.02(3) 


2.0 


2.0 


2.0 


0(4) 


B 


1.6(3) 


2.4(3)" 


1.35(5) 


1.73(9) 




n 


0.96(2) 


1.0 


0.92(2) 


1.0 


0(5) 


B 

n 


1.6 

0.04(1) 










a 


3.8172(1) 


3.8206(1) 


3.8231(1) 


3.8198(1) 




b 


3.8822(1) 


3.8851(1) 


3.8863(1) 


3.8849(1) 




c 


11.6707(4) 


11.6757(4) 


11.6809(2) 


11.6762(3) 



" Equivalent isotropic B from anisotropic results. 

YBajCUjO, - Orthorhombic model - no opproximotions 



{Cu(l)}0(l)A4). 
Cu(l) {Cu(2)}0(l)0(2)A3)2 

{Ba}0(l)A2)A3)A4). 
Cu(2) mO(2)A3l 

{0(l)}Cu(l)Cu(2)Ba, 

{0(2)} Cu(2),Ba,Y, 




Cu(2) 



{0(3)} Cu(2),Ba,Y, 
{0(4)}Cu(l),Ba, 



YBOjCUsOj - Tetragonal 
Cu(l)Cu(2)2BajYO(l),0(23)4 




CU(1) 



Cu(2) 



Cu(2] 



{Cu(l)}0(l), 
{Cu(2)} 0(1)0(23), 
{Ba}0(l)A23), 
{Y} 0(23)e 

{0(l)}Cu(l)Cu(2)Ba, 
{0(23)} Cu(2)2Ba2Y2 



Fig. 20. Formula units, bonding schemes and atomic coordinations in 
the structures of orthorhombic and tetragonal BaaYCujOT and 
BajYCujO,,. 



(2.33 +) for X = 7.0. Tliis means that adding oxygen to 
tlie 0(4) sites on tlie plane of the Cu(l) atoms increases 
not only the oxidation state of Cu(l), but also that of 
Cu(2), i.e., a charge transfer between the two copper 
atoms takes place during the oxidation, or reduction, of 
the compound. The Ba-0 (1) and Cu(2)-0(23) bonds at 
X = 6.0 are strained because they are incommensurate 
when the oxidation states of Cu(l) and Cu(2) are 1+ and 
2+. As the oxidation progresses, however, the incom- 
mensurability of these bonds decreases, and the strains 
are almost completely relieved atx = 7.0. This behavior, 
therefore, suggests that the observed charge distribution 
(and all the properties related to it) at different oxygen 
stoichiometrics is a consequence of the commensurabil- 
ity of the Cu-0 and Ba-0 bonds. 

3.2.5 The System RBajFcjOs+x 

The substitution of copper by other metal species in 
the structure of YBaaCujG ^ has been used extensively to 
investigate the correlation between the changes induced 
into the structure by the substituting cations and the 
variation of the superconducting properties. The first 
studies in this direction involved the partial replacement 
of either the chain copper atoms Cu(l), e.g., [61], or the 
plane copper atoms Cu(2), e.g., [62]. From these exper- 
iments it was found that the introduction of dopants 
decreases the value of Tc and that the change is very 
drastic when the atoms Cu(2) are replaced even by small 
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amounts of other elements such as Zn and Ni. These 
results were taken as evidence that the two-dimensional 
Cu(2)0(23)2 layers play a fundamental role in the super- 
conducting process. Full replacement of copper by other 
elements, therefore, is important for understanding why 
copper, of all transition metals, seems to be necessary 
for superconductivity to occur in this structural type. 
For this reason collaboration between the NBS reactor 
and the University of Oslo was established in the early 
nineties to study the system in which copper is entirely 
replaced by iron. 

The first compound investigated in this series of 
experiments was stoichiometric YBa2Fe308 [63]. The 
nuclear structure of this material has the symmetry of 
space group PAImmn and is illustrated in Fig. 21. The 
configuration of the atoms in the unit cell is very similar 
to that of the superconductor YBa2Fe307, with the 
exception that the iron atoms Fe(l), corresponding to 
Cu(l), have octahedral coordination, rather than four- 
fold planar, and the Fe(l)06 octahedra are arranged in 
layers. This significant difference is a consequence of 
the fact that all possible oxygen sites on the Fe(l) plane 
are fully occupied, giving an oxygen content of eight 
atoms per formula unit. A second consequence is 
that the Ba atoms have twelve-fold cuboctahedral coor- 
dination. This configuration results in a rather strained 
structure in which the polyhedron around Ba is 
considerably compressed. 

During the structural refinements based on neutron 
diffraction data, it was realized that the nuclear 
structure of the iron compound did not account for all 
the observed intensities, and especially for those of 
some strong reflections at low values of the 2 6 angle. 
Polarized neutron diffraction experiments revealed that 
these intensities have magnetic origin, and a subsequent 
analysis showed that they are consistent with a mag- 
netic structure in which the iron moments lie on the a, 
b plane and are coupled antiferromagnetically within 
each Fe02 layer, as well as along the c-axis, as 
indicated by the arrows associated with the Fe atoms in 
Fig. 21. The magnetic moments of the two iron atoms 
were found to be equal to each other (3.49 //-b in both 
cases) thus showing that the oxidation states of Fe(l) 
and Fe(2) are also the same (3h-). The results of our 
study were confirmed later by x-ray powder diffraction, 
magnetic susceptibility and Mossbauer spectroscopy 
[64, 65]. In another experiment, however, the differ- 
ences between the magnetic hyperfine fields at Fe(l) 
and Fe(2) were interpreted as evidence that the valence 
states and the magnetic moments of the two iron 
sites are not the same [66]. In order to resolve this 
uncertainty we have carried out a complete bond- 
valence analysis of the structure and found that the 
atomic configuration of the iron compound is consistent 
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Fig. 21. View of the tetragonal crystal structure of YBa2Fe30g. To 
facilitate comparison, the labeling of the atoms is consistent with that 
of Fig. 19. The directions of the magnetic moments of the two crystal- 
lographically independent iron atoms are indicated by arrows. 



with an equal distribution of charge on the Fe(l) and 
Fe(2) sites, thus confirming our original conclusion 
based on the determination of the magnetic structure 
[67]. 

Soon after the original study described in Ref. [63], 
it was realized that YBa2Fe308+x can exit in a narrow 
range of compositions, with -0.2 <x< 0.1. The location 
of the oxygen in excess of the Og stoichiometry was 
determined in the analysis of the Ca-doped compound 
(Yi_, Ca,,) YBa2Fe308« [68]. As indicated in Fig. 21 
these atoms occupy the 0(int.) sites, on the plane of the 
Y atoms. Their presence in the structure introduces 
significant disordering of the atoms in the 0(23) sites 
and increases the average coordination of Y and Fe(2). 
On the other hand, the presence of oxygen vacancies 
involves the 0(45) sites and cause a decrease of the 
average coordination of Ba and Fe(2). 
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3.2.6 Hg-Based Superconductors 

In the early nineties, superconductivity with 
rc=94K was discovered in the mercury based com- 
pound HgBa2Cu04+t [69]. The structure of this material 
was first determined by x-ray powder diffraction and is 
schematically illustrated in Fig. 22. The atomic con- 
figuration of this compound can be conveniently 
described by means of the layer sequence 

. . . [(BaO),(HgO,)o(BaO),(Cu02)o](BaO), . . . 

where the square brackets include the content of one 
unit cell and the subscripts c and o indicate if the cation 
is at the center or at the origin, respectively, of the mesh 
of each layer. The formulae in parentheses indicate the 
chemical composition of the layers. Since the extra 
oxygen atoms 0(3) located on the plane of the Hg atoms 
are necessary to increase the oxidation state of copper, 
it is important to know how many of them can be in- 
corporated in the structure and where exactly they are 
located on the layer. As we have mentioned earlier, 
X rays are not sufficiently sensitive to small amounts of 
oxygen in the presence of heavy atoms such as Ba and 
Hg, and for this reason the compound was studied at the 
NIST reactor in a collaborative project with the French 
and Russian authors who discovered this important 
material [70]. 

The copper atoms in this structure (see Fig. 22) have 
a bi-pyramidal coordination similar to that found for 
Lai85Sro.i5Cu04, with an apical Cu-0(2) distance of 
2.78 A, compared with the in-plane distance Cu-O(l) of 
1.94 A. Most of the Ba atoms have eight-fold coordina- 
tion. However, those located below and above the atoms 
0(3) are nine-coordinated and the coordination poly- 
hedron may be considered a mono-capped square 
antiprism. Similarly, the Hg atoms are generally two- 
coordinated, except in those unit cells in which the sites 
0(3) are occupied by oxygen atoms. 

The compound HgBa2Cu04+;t is the first member of a 
homologous series of formula HgBa2R„_iCu„02„+2+t and 
having a structure represented by the layer sequence 

. . . [(BaO),(Hga)„(BaO),(Cu02)o{(«-l)(R)c(CuO,)„|](BaO), . . . 

where R is an element situated on a layer that does not 
contain oxygen atoms and is therefore similar to the 
yttrium layer in YBa2RCu30t. Thus, the layer sequence 
of the compound HgBa2Cu206-x, corresponding to 

n = 2 is 

. . . [(BaO),(HgO,)o(BaO)e(Cu02)o(R)c(Cu02)o](BaO)e . . . 

and so on. The sequences show that the struc- 
tures of these compounds are made by blocks 




Fig. 22. Unit cell of the crystal structure of the 94 K superconductor 
Hg Ba2Cu04+,. In this cell, oxygen sites 0(3) are only partially occu- 
pied. The copper atoms are in a bi-pyramidal coordination similar to 
that of Fig. 17. The Ba atoms are in eight-fold coordination in those 
unit cells in which 0(3) is empty and in nine-fold coordination in the 
others. Similarly, Hg may have two-, three-, four-, five-, and six-fold 
coordination, depending on how the 0(3) sites occupied by oxygen 
are clustered together. 



. . . (Cu02)o(R)c(Cu02)o . . ■ having the configuration of 
perovskite, and by blocks . . . (BaO)c(HgOt)o(BaO)c . . . 
with the rock-salt structure. To avoid writing complex 
formulae, these materials are usually identified by 
indicating the number of the Hg, Ba, R, and Cu cations, 
respectively. Thus, the first member of the series is 
1201, the second 1212, etc.. All members have the 
tetragonal symmetry of space group PAImmm 
and approximate lattice parameters a ~ 3.85 A and 
c ~ 9.5 -I- 3.2 (n-1) A. The oxygen content x depends 
on the preparation method and the annealing conditions 
of the sample. The maximum value of x in each member 
of the series, however, seems also to be function of n , 
i.e., of the number of layers that constitute the 
perovskite block [71, 72]. The critical temperature T^ is 
strongly dependent on x . An example of the behavior of 
Tc is provided by the 1201 compound and is illustrated 
by the graph of Fig. 23, taken from Ref. [73]. 
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Fig. 23. Plot of Tc versus the composition x for the superconductor 
Hg Ba2Cu04„. 



For X < 0.06, Tc is equal to zero, then it increases with 
increasing x, reaches a maximum of Tc = 94 K for 
X ~ 0.17, and then decreases. Similar behavior has been 
observed also for other members of the series. For each 
of these compounds, in which R = Ca, the maximum 
value of Tc is a function of n , and, more specifically, 
it increases for increasing n up to n = 3, and then 
decreases as shown Table 3. 



Table 3. Composition and T^ values in the homologous series 
HgBa2Ca„_,0„+2+, 



Compound 




X 


T, 


HgBa2Ca5CU(,0,4« 


1256 


0.065 


126-128 K 


HgBajCuGi^ 


1201 


0.065 


94-95 K 


HgBa2CaCu20(,« 


1212 


0.22 


126-128 K 


HgBa2Ca2Cu308„ 


1223 


0.41 


133 K 


HgBa2Ca5CU(,0,4« 


1256 


0.40 


114K 



4. Zeolitic Materials 

(This section was written by B. H. Toby, Leader, 
Crystallography Team) 

4.1 Introduction 

Zeolites and related microporous materials (herein 
referred to as zeolitic materials) are crystalline materials 
consisting of alternating regions of framework atoms 
and voids. The framework is typically composed of 
tetrahedral aluminum and silicon atoms linked by 
oxygen atoms and the voids or pores can range to just 
accommodating a water molecule, to pores large enough 
to fit a Ceo molecule or DNA fragment. 



Zeolitic materials frequently form highly symmetric 
and complex structures that are elegant to behold. 
Further, these materials are industrially important for a 
number of diverse applications. At least one zeolite- 
based catalyst, and likely more than one, is employed to 
produce all domestic gasoline. Zeolitic catalysts are in 
creasingly used for the preparation of pharmaceuticals 
and other fine chemicals. Zeolitic materials have novel 
capabilities, which can often be enhanced to perform a 
particular chemical separation. They are frequently 
used commercially to separate gases, for example 
producing N2 or O2 from air. Another application is 
to absorb water, for example to keep double-pane 
windows from fogging. 

Most zeolitic materials have negatively charged 
frameworks that require cations for charge balance. 
This allows these materials to be used for ion exchange. 
In our homes, zeolites are formulated in detergents to 
remove the calcium ions that make water hard. Differ- 
ent materials are used to sequester radioactive ions for 
environmental remediation. 

Zeolitic materials rarely crystallize as single crystals. 
This means that most structural studies must utilize 
powder diffraction. Further, these materials frequently 
contain both light and heavy atoms and thus neutron 
scattering information is essential to determine sitting 
of light atoms. The use of both neutron and synchrotron 
x-ray diffraction data in many cases has been a pre- 
requisite for accurate structural determinations (for 
example, see [74]). 

The following paragraphs summarize results from 
many recent studies performed at NIST where neutron 
diffraction played a pivotal role. 

The 32 detector BT-1 neutron powder diffractometer 
in the NIST Center for Neutron Research allows the 
instrumental resolution to be tailored to suit the needs of 
the sample[75]. This has been very effective for the 
study of zeolitic materials, since optimum resolution is 
needed at relatively low Q values, as opposed to dense- 
phase materials. The instrument was further improved 
for zeolitic materials in the mid-1990s, when the 1.54 A 
Cu(220) 75° takeoff angle monochromator was 
replaced with a Ge(311) monochromator with a wave- 
length of 2.078 A. Data can be collected approximately 
twice as fast with Ge(311) compared with the Cu(220) 
monochromator and resolution is even better tailored to 
the needs of zeolitic diffraction experiments. 

4.2 Zeolite RHO 

The zeolite RHO framework is shown in schemati- 
cally in Fig. 24. This material has been of interest in part 
due to the extreme flexibility of the framework, as 
shown in Fig. 25. Cesium-exchanged zeolite RHO was 
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S8R 



Fig. 24. Schematic representation of the Zeolite RHO framework. 
Vertices indicate tetrahedral atoms (oxygen atoms are omitted for 
clarity). The single eight-ring (S&R) pore opening is indicated. 
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Fig. 25. Configuration of the Single Eight Ring in Zeolite RHO 
materials. Tetrahedral atoms are shown as open circles; oxygen atoms 
are shown as larger filled circles. The A parameter indicates the 
degree of distortion. 



one of the first materials to be studied in the early 1970s 
on the then new five-detector BT-1 instrument [76]. 
The framework topology has continued to be of 
interest. In one recent study, an aluminosilicate RHO 
material, Li76Csi3Na2oAlii4Si366096, was com- 
pared to a novel aluminogermanate RHO material. 



Lii3.9Cs5 24Nao.24Al24Ge24096 [77]. This study demon- 
strated that despite having the same framework, the 
different compositions resulted in different cation sites. 
In these systems, Li sites could not be determined 
without use of neutron diffraction. 

The high degree of flexibility in the RHO framework 
leads to what is known as the trap door effect with some 
cations, such as Cd. At room temperature, the cation 
binds in the single-eight ring window (S8R, see in [24]) 
but at elevated temperatures, the cations migrate to 
higher coordination sites. Initial reports indicated that 
when the material was cooled, the cations returned to 
the S8R sites [78, 79]. When the cations are located in 
the S8R site, they block access to the zeolite pores. Thus 
it is possible to trap species in this material, by sorbing 
the guest at elevated temperatures. When the tempera- 
ture is lowered, the guest is trapped when the cations 
return to the S8R sites. Subsequent study indicated that 
cations are stable in the S8R site only when hydrated 
and that H2O (or possibly OH ) remains bound to the 
cations at much higher temperatures that previously 
understood [80]. With heating, the water desorbs and 
the cation migrates to a higher coordination site. The 
cation remains in this higher coordination site when the 
temperature is reduced, except in the presence of water 
vapor. Thus, the trap door effect is actually a hydration/ 
dehydration phenomenon rather than simply a tempera- 
ture-mediated effect. These results also call into ques- 
tion the observation of negative thermal expansion in Sr 
exchanged RHO [81, 82]. It is now believed the con- 
traction of the unit cell upon heating is due to a re- 
versible loss of water. 

4.3 Structure-Directing Agents in Synthesis of 
CON Frameworli Materials 

Many zeolitic materials are prepared from gels con- 
taining a sacrificial organic amine, called a structure-di- 
recting agent (SDA), that is trapped in the pores during 
synthesis. To free the pore spaces, the SDA is oxidized 
away by calcinations in air. The influence of the SDA 
species during crystallization is poorly understood. In- 
deed a single SDA can be used to produce different 
materials from the same gel composition, by variation 
of processing conditions [83]. A study combining neu- 
tron and synchrotron diffraction with molecular me- 
chanics was recently undertaken to investigate how 
SDA cations interact with a set of zeolitic materials 
[84]. The borosilicate zeolitic materials, SSZ-33 and 
CIT-1 have related structures and have been designated 
the structure code CON by the Structure Commission of 
the International Zeolite Association. Unlike CIT-1, 
which is nearly free of defects, SSZ-33 has a very high 
density of stacking faults (>30 %) [85]. Several organic 
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cations, such as species SDA 1 shown in Fig. 26, can be 
used to synthesize SSZ-33. However, only SDA 2 has 
been found to synthesize CIT-1. Interestingly, cation 
SDA 3, while quite similar to SDA 2, cannot be used to 
prepare any material in the CON family. 




SDA 1 



(CH3)3N^ 



(CH3)3N: 




SDA 2 



(CH3)3N^ 




SDA 3 



Fig. 26. Structure-directing agents studied for preparation of ^ON 
framework materials. 



The location of the SDA cation in CIT-1 prior to 
calcinations was determined using a simultaneous 
Rietveld fit to both neutron powder diffraction data and 
synchrotron x-ray powder diffraction data. Determi- 
nation of the cation orientation was made possible by 
selective deuteration of the methyl groups of the quater- 
nary nitrogen. Since crystallographic methods do not 
prove uniqueness of the resulting model, molecular 
modeling was used to tabulate all possible sites where 
the SDA cation could be accommodated in the CIT-1 
pores. Based on these computations, it was concluded 
that the crystallographic model is the only plausible 
result. This model, shown in Fig. 27, demonstrates the 
very tight fit between the SDA ions and with the CIT-1 
framework. 

Subsequent molecular modeling work was performed 
using the crystallographic model as a starting 
point. These results confirmed that four molecules 
of either SDA 1 or SDA 3 can pack in the CIT-1 
pores without any unfavorable energetic inter- 
actions and with similar energetic. However, packing 
SDA 3 in the same voids required significantly more 
energy, due to Van der Waals repulsions. This suggests 




Fig. 27. Packing of SDA-2 in CIT- 1 . Carbon atoms are shaded and 
hydrogen atoms are shown in white, both at Van der Waals radii; 
nitrogen atoms cannot be seen. The CIT-1 framework is shown as 
dark lines, where oxygen atoms have been omitted. 



why SDA 3 does not synthesize the CON framework, 
since it cannot pack effectively to form the same voids. 
Molecular modeling of SDA packing in the presence of 
stacking faults showed no significant energetic differ- 
ences between SDA 1 and SDA 2, indicating that these 
stacking faults likely arise as a kinetic effect rather than 
due to thermodynamics. 

4.4 Novel Lithosilicate Molecular Sieves 

There is great interest in finding new zeolitic frame- 
works, e.g., materials with new pore structures, 
since these materials have the potential to offer unique 
properties. For this reason there was considerable 
excitement accompanying the discovery that novel 
microporous materials can be prepared with frame- 
works composed of Si04 and Li04 tetrahedral [86, 87]. 
The Li04 tetrahedral appear to be more flexible than 
Si04 or AIO4 tetrahedral and thus these lithosilicate 
materials contain building units that would be highly 
strained in an aluminosilicate. The framework of one 
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such material, RUB-29 is shown in Fig. 28. The approx- 
imate framework structure was determined using single- 
crystal synchrotron diffraction from an extremely small 
(—200 jjim^) crystal. Neutron diffraction was then used 
to determine the Li sitting. NMR measurements indicate 
that RUB-29 is also novel in that all Li atoms in this 
material appear to be mobile below 250 °C. This 
suggests that lithosilicates could be useful as ion 
conductors, perhaps in batteries or fuel cells. 




Fig. 28. Framework geometry of litho silicate RUB-29. Oxygen 
atoms are shown as red spheres. The Si04 units are shown as brown 
tetrahedra and Li04 units as pink tetrahedra. Charge balancing 
cations, sited in the pores, are omitted for clarity. 



5. Crystallography in Engineering 
Research 

(This section was written by Thomas Gnaupel-Herold, 
NIST Center for Neutron Research) 

5.1 Introduction 

An important role of mechanical engineering is the 
calculation of safety margins and the lifetime prediction 



of parts and components. A key factor for these consid- 
erations is the knowledge of stresses a component would 
experience in service. While the fraction of external 
stresses is usually well-known, residual stresses can be 
problematic because they superimpose on the external 
stresses resulting in a total stress level that can exceed 
the strength of the material, thus causing unexpected 
failures. In addition to manufacturing induced residual 
stresses, there are also service-induced stresses originat- 
ing from fatigue, corrosion and other processes that 
often determine the lifetime of an engineering compo- 
nent. 

Due to the high penetration of neutrons in many 
materials, neutron diffraction residual stress analysis 
(NDRSA) is one of the few methods able to provide 
information about the three-dimensional strain and 
stress distribution at depths of several centimeters in 
engineering components nondestructively. In compari- 
son, X rays from laboratory sources in the energy range 
5 keV to 20 keV are limited to about 20 |xm penetra- 
tion. Thus, NDRSA can act as a tool for providing basic 
data for the design of a component, model validation of 
finite element calculations or even for quality control at 
certain stages of the life of a part. Recognizing these 
uses and the need for such a tool for American industry 
and academia, a collaborative program of research in 
this area was begun in the early 1980s. The NIST Center 
for Neutron Research began construction of a new, 
dedicated double axis system for residual stress, texture 
and single crystal analysis (DARTS) in 1991. The instru- 
ment has been in operation successfully since fall 1995. 

5.2 Principle of the Method 

As diffraction is basically a length measurement, lat- 
tice spacing act as a strain gauge for a stress in a certain 
direction of the part or component. The relationship 
between the wavelength of the neutrons A, the lattice 
spacing di,ti of the atomic planes described by the Miller 
indices hkl, and the diffraction angle is given by Bragg's 
law. 



A = 2diM sin 



(13) 



The fif- spacing information does not come from a point, 
it describes the volume average over a region called the 
sampling volume. The shape and dimensions of the 
sampling volume are prescribed by two slit systems. 
The first set of slits defines the incident beam which is 
diffracted throughout its entire path through the sample. 
The secondary slits in front of a neutron detector single 
out a narrow part of that diffracted beam. The intersec- 
tion of both beams defines the sampling volume (see 
Fig. 29). Usually the wavelength is chosen so that 26hki 
is approximately 90° because then the cross-section of 
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Fig. 29. Schematic of diffraction geometry and slit arrangement. 



the sampling volume becomes rectangular. Then, its 
dimensions are the width of primary beam times the 
height times the horizontal opening of the receiving slit. 
The position of the sampling volume is constant. Three- 
dimensional mapping of strain in a component is ac- 
complished by using an x-y-z- positioning table which 
can move any point in the component into the sampling 
volume. 

If a reference value for the cf -spacing can be found, 
e.g., by measuring li-spacings in a small coupon in 
which all long-range stresses are relieved, so that 
the calculation of strain from the measured cf-spacings 
becomes straightforward 



d.. 



s(hkl, ip, ili)= 

Uhkl 



(14) 



The angles cp and ip describe the orientation of the 
component (specimen) with respect to the reference 
frame. Using strain measurements in at least six 
independent directions, the complete strain tensor can 
be obtained from 



e(hkl, (p, ip) = ZjBi 



niti. 



(15) 



in which n = (cos ^ sin t//, simp smij/, cosi//). In order to 
obtain stress from strain, the elastic moduli in Hooke's 
law have to replaced by the plane specific values Ehu 
and Vkki for translating lattice strain into macroscopic 
stress. 



5.3 Instrumentation 

A residual stress neutron diffractometer is nothing 
more than a specialized powder diffractometer. The 
differences are in the sophisticated slit assembly, the 
available wavelength range, the resolution and, most of 
all, an x-y-z -table for moving and positioning heavy 
specimens. 

A residual stress measurement requires slit systems 
that define the sampling volume because this is the best 
way to assign the measured (i-spacing to a point or, 
because it is not possible to measure a point, to a small 
region within the sample. Without these slits, the 
d-spacing would represent the spatial average over the 
entire illuminated region. Also, because of the large 
penetration of neutrons, this region would extend over 
several cm"* instead of a few mm^ as required by most 
problems. The available neutron flux usually makes 
sampling volumes <0.5 mm^^ unfeasible. 

Because of the requirement of having 2 dhu around 90° 
for good spatial resolution, there is no need for a large 
angular opening in the scattering plane. At NIST, the 
detector is a linear position sensitive detector of 100 mm 
length which under normal conditions provides an 
opening of about 8° — enough to record a complete 
diffraction peak. 

In order to bring almost any desired peak to 90° the 
instrument has a total of three different monochroma- 
tors with six available reflections. Together with a vari- 
able monochromator diffraction angle Bu, the wave- 
length range is from 0.95 A to 3.13 A. There is also 
enough overlap between the wavelength range of each 
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monochromator so that the optimal setting can be 
chosen that is given by the figure of merit (FM). The 
FM is, roughly speaking, the neutron flux at the sample 
position divided by the square of the standard deviation 
of the recorded diffraction peak. It is a measure of the 
counting time necessary to reach a certain for the 
chosen d-spacing. In fact, a residual stress machine is 
optimized for the best possible FM which makes it more 
of a medium resolution diffractometer. Unlike powder 
diffractometers which may be more tuned towards 
better resolution and make up for the loss of intensity by 
bigger samples, bigger beams, and multiple detectors, 
the design of a residual stress diffractometer is driven by 
the size of the sampling volume. The goal is to measure 
a strain with a given accuracy for a volume ~ 30 mm' or 
smaller in the shortest possible time. 

Big improvements in the figure of merit were made 
possible by the recent development of perfect silicon 
double focusing monochromators. The concept is to 
utilize the large in-pile acceptance (area about 
18cm^Xl3cm^) and to focus otherwise unused neu- 
trons into a narrower spot at the sample position. The 
last generation of these devices utilizes packets of thin 
silicon wafer blades that are stacked vertically and tilted 
slightly with respect to one another so that the whole 
assembly has a fixed vertical curvature. Horizontally, 
i.e. within the scattering plane, the curvature is 
produced by elastic bending of the packets in a four- 
point bending device. Figure 30 shows a schematic 
diagram of the device used at NIST. 



The optimal horizontal curvature, i.e., the curvature 
that produces the highest FM, depends on a variety of 
factors but the most important one is the monochroma- 
tor diffraction angle 0m- Because the versatility of a 
variable du is needed for changing the neutron wave- 
length it is also necessary to change the horizontal 
curvature to the respective optimal setting. Such a 
device is used at the DARTS diffractometer at the 
NCNR. 

Aside from the instrumentation for delivery and 
detection of neutrons, there is also peripheral equip- 
ment. This includes a four-circle goniometer for single 
crystal and preferred orientation measurement, a load 
frame for the investigation of elastic and plastic 
deformation characteristics and heating devices for the 
measurement of stress at elevated temperatures. 

5.4 Examples of Neutron Diffraction Residual 
Stress Measurements 

Most of the residual stress measurements are related 
in one way or another to industrial applications. There- 
fore the most frequently investigated samples are the 
technologically important materials such as steel, 
aluminum, nickel-alloys and ceramics. The diagram of 
Fig. 31 shows an overview of the materials that were 
investigated 1999 at the residual stress diffractometers 
at NIST, HIFR (Oak Ridge), ILL (France) and ISIS 
(UK). 




Fig. 30. Schematic drawing of the double focusing perfect silicon monochromator. On the right side, the useful 
reflections of the device are shown. 
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Fig. 31. Materials investigated with residual stress diffractometers world wide (selected). 



NDRSA provides depth profiling capabilities at 
spatial resolutions that match both the requirements 
imposed by finite element meshing and by the 
characteristic length scale within which steep strain 
gradients are present in a sample. This way even detailed 
maps of three dimensional strain and stress distributions 
can be obtained. 

Examples for the work that has been done at the 
NIST Center for Neutron Research include welds in 
aluminum, steel and stainless steel; automotive springs; 
ceramics and metal-matrix composites; induction- 
hardened fatigue specimens; structural parts of aircrafts. 
A good example for the strain mapping capabilities of 
NDRSA is a project done in collaboration with the U.S. 
Department of Transportation and residual stresses in 
rails. The purpose of this investigation was to determine 
the effect of grinding strategies on the distribution and 
magnitude of service-induced residual stresses. Their 
importance relies on the fact that they are critical 
parameters in the estimation of the growth rate of 
fatigue defects. This information is used to assess rail 
inspection frequencies to assure defect detection prior 
to catastrophic rail failure. Figure 32 shows the two 
dimensional residual stress distribution a unworn and a 
worn rail obtained from a mesh of sampling volumes of 
3 mm^ X 3 mm^ X 3 mm^ (right side). 



The above example illustrates current capabilities of 
NDRSA. Stresses can be obtained from depths of about 
3 cm in steel, 6 cm in aluminum and 2 cm in nickel 
using sampling volumes of 30 mm^ to 50 mm^ Further 
instrumental improvements will not increase that depth 
significantly because of the exponential attenuation of 
the neutron beam. The gains in the figure of merit will 
rather go into improved stress accuracy and smaller 
sampling volumes. For example, an improvement in the 
figure of merit by a factor of two enables a sampling 
volume half the previous size at the same counting time 
and stress accuracy. 

6, Search for Novel Superconducting 
Materials 

(This section was written by Q. Huang, NIST Center for 
Neutron Research) 

As mentioned above, a lot of work has been done on 
the structural determinations of the high T^ supercon- 
ducting materials using neutron powder diffraction. In 
particular, materials with the oxide perovskite type 
structure have attracted interest in recent years because 
this structure type seems to be an excellent structural 
framework for promoting superconductivity, eclipsing 
the intermetallic compounds which have been the source 
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Fig. 32. Residual stress distribution in an unknown (top) and a worn (bottom) rail. Red colors indicate tensile stresses, blue stands 
for compressive stress. Values on the color bar are given in megapascals. The right side shows the mesh of sampling volumes with 
size and location within the rail. 



of many superconducting materials in the past. How- 
ever, the recent discovery of superconductivity in MgB2 
[88] suggests that intermetallic compounds with simple 
structure types are worth serious reconsideration as 
sources of new superconducting materials. This sugges- 
tion is corroborated by the discovery of superconductiv- 
ity at 8 K in the intermetallic compound MgCNi3, 
which has the perovskite structure [9]. The high propor- 
tion of Ni in the compound suggests that magnetic 
interactions may be important to the existence of the 
superconductivity. MgCNis, in fact is best considered 
as the three-dimensional analog of the layered nickel 
borocarbides, typified by LuNi2B2C, which has a T^ 
near 16 K [89]. The discovery also indicates that MgB2 
and MgCNis will not be the only superconductors of this 
kind, and that a significant new class of superconducting 
materials may be found down this path in the future. 
Figure 33 shows the structure of the non-oxide 
perovskite intermetallic compound MgC^Nis which 
has a Pm3m symmetry with atomic positions of 
Mg: (0, 0, 0), C: (1/2, 1/2, 1/2), and Ni: (0, 1/2, 1/2). 
The C atom is surrounded by Ni atoms forming CNie 
octahedra which, in the oxide perovskites, are formed 
by a metallic atoms with six surrounding oxygen atoms. 
The magnetization data [89] show that the magnetic 
onset for the superconducting transition ranges between 
7.1 K for a nominal carbon content of 1.1 per 
MgCXNi3, to 7.4 K for nominal carbon content x= 1.5. 
For nominal C contents between x = 1.1 and x = 1.0 the 
superconducting transition turns off abruptly. The C 



content x in the structure must play a major role for the 
superconductivity to exist. 




iMg ONi •C 



Fig. 33. Schematic representation of the crystal structure of the 9 K 
superconductor MgCNis. 



The structure refinement, using the neutron powder 
diffraction data collected at BTl, reveals that the com- 
position of the superconducting phase for nominal 
x= 1 .25 (T, = 13 K) is MgC(,.96Ni3. This result is in 
agreement with the small amount of unreacted graphite 
(2 % by weight) found in the sample. Therefore the data 
suggest that the highest T^ in the perovskite is obtained 
for the stoichiometric composition MgdoNij. Further 
the data suggest that the superconductivity disappears 
abruptly for carbon contents between 0.96 and approxi- 
mately 0.90 per formula unit. 
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Fig. 34. Plot of the a -parameter (top) and the thermal factors (bottom) as function of temperature, in MgCNis. The behavior 
of these parameters shows no detectable anomalies at the critical temperature. 



MgCo.96Ni3 shows no magnetic or structural transi- 
tions between 2 K and 295 K [90] . Fitting with a 
Gausian function all the peaks above 70° 2-theta at 
different temperatures indicates that the peak width and 
shape are close to the instrumental resolution, i.e., no 
observable particle and/or strain broadening are present 
in the measured temperature range, and that no evidence 
of other structural distortions was found with decreasing 
temperature. The refined lattice parameter a and Debye- 
Waller factors decrease smoothly as the temperature 
decreases. Careful measurements near T^ (7.3 K) show, 
as shown in Fig. 34, that the changes of those para- 



meters are within the experimental uncertainties, which 
are relatively small. Therefore if there is any structural 
anomaly at T^, it would be undetectable with the degree 
of precision of these measurements, which is of the 
order of 0.01 %. 
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